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Zastosowanie metod programowania genetycznego 
w procesie maksymalizacji wydobycia węglowodorów 
przy zastosowaniu symulatora złożowego 


Artykuł poświęcono zastosowaniu metody programowania genetycznego dla celów optymalizacji wydobycia ropy nafto- 
wej na przykładzie testowego złoża węglowodorowego. Prezentowane zagadnienie optymalizacyjne jest prostym przykła- 
dem problemu optymalnej kontroli i polega na doborze wydajności wydobycia ropy naftowej w przyjętych przedziałach 
czasowych w taki sposób, aby w zadanym całkowitym czasie eksploatacji uzyskać maksymalne wydobycie sumaryczne 
przy minimalnym wydobyciu wody. Problem rozwiązano przy zastosowaniu algorytmu genetycznego, kodującego dozwo- 
lone wartości wydajności wydobycia z listy wartości dozwolonych. Z jednej strony działanie takie jest charakterystyczne 
dla metod programowania genetycznego, zaś z drugiej redukuje istotnie przestrzeń rozwiązań. W artykule zastosowano 
algorytm genetyczny Hollanda, dla którego zaimplementowano krzyżowanie wielopunktowe oraz adaptację prawdopodo- 
bieństw krzyżowania i mutacji na podstawie tzw. współczynnika zróżnicowania populacji. Działanie tak zdefiniowanego 
mechanizmu adaptacji jest następujące: jeżeli zróżnicowanie populacji rośnie, liniowo zwiększane jest prawdopodobień- 
stwo krzyżowania, a zmniejszane prawdopodobieństwo mutacji; w przeciwnym wypadku (zróżnicowanie populacji male- 
je) działa mechanizm odwrotny, tzn. zmniejsza się prawdopodobieństwo krzyżowania, a zwiększa prawdopodobieństwo 
mutacji. Taka metoda z jednej strony gwarantuje różnorodność populacji, z drugiej zaś zapewnia dobrą eksploatację prze- 
strzeni rozwiązań. Przeprowadzono szereg testów mających na celu zweryfikowanie efektywności algorytmu w zależno- 
sci od liczby punktów krzyżowania (krzyżowanie 1-, 2-, 3-punktowe) oraz długości chromosomu. Wykonane testy wska- 
zują na zadowalającą zbieżność algorytmu, niezależnie od wartości badanych parametrów. Przyjęcie funkcji w określonej 
postaci spowodowało premiowanie przez algorytm niższych wartości wydobycia, co wynika z nieliniowego przyrostu wy- 
dobycia wody dla wyższych wartości wydobycia ropy naftowej. 


Słowa kluczowe: algorytmy genetyczne, programowanie genetyczne, optymalizacja wydobycia, symulacje złożowe. 


Application of Genetic Programming Methods for the Optimization of Hydrocarbon Production 
by using a Reservoir Simulator 


The paper addresses the problem of oil production optimization by genetic programming methods. The specific example 
of the problem presented in the paper belongs to the class of, so called, optimal control problems. It consists in finding the 
time variable rates of oil production that result in the maximum of the total oil production while keeping the total water 
production at a minimum available level. The problem is solved by a genetic algorithm, that assumes the production rates 
from the list of the allowable values. This approach typical for genetic programming methods significantly reduces the space 
of possible solutions. The article uses the Holland genetic algorithm for which multi-point crossing has been implemented 
and the adaptation of crossing and mutation probabilities based on so the called coefficient of population variability. The 
adaptive mechanism makes the crossing probability increase and mutation probability decrease for population variability 
increasing with time, while the crossing probability decrease and mutation probability increase for the variability decreas- 
ing with time. This mechanism guarantees the population variability to be at on appropriate level and at the same time, the 
extrapolation process for the solution space to be effective. Several tests were performed to verify the actual effectiveness 
of the algorithm for various number of crossing points (1, 2, 3 — crossing points) and chromosome length. Their results 
show a satisfactory convergence of the method to the final solution independent of the varying parameters values. Adopting 
a function in a specific form resulted in an algorithm for lower mining values, resulting from a nonlinear increase in water 
extraction for higher oil production values. 
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Istnieje wiele problemów, których rozwiązanie klasycz- 
nymi metodami skazane jest na niepowodzenie. Niektóre 
z nich powstały jako zagadnienia teoretyczne, jednak ich 
efektywne rozwiązania posiadają istotne znaczenie prak- 
tyczne w wielu dziedzinach nauki i techniki. Przykładem 
może być problem „komiwojażera” czy „plecakowy”. Nie- 
doskonałość klasycznych metod optymalizacji w zastoso- 
waniu do przedstawionych zagadnień spowodowała, że 
rozpoczęto poszukiwania metod działających według in- 
nych zasad. Szukając inspiracji, skierowano uwagę na świat 
przyrody i okazało się, że ewolucja rozwiązała już takie 
problemy. Jednym z najbardziej popularnych zastosowań 
rozwiązań ze świata biologii wykorzystanym przy wyja- 
śnianiu zagadnień optymalizacyjnych są algorytmy gene- 
tyczne [18] mające źródło w teorii ewolucji gatunków i za- 
liczane dzisiaj do rodziny ewolucyjnych metod optymali- 
zacyjnych [3, 24, 25]. Spośród wielu zastosowań algoryt- 
mów genetycznych na zainteresowanie zasługują metody 
tzw. programowania genetycznego spopularyzowane przez 
Johna Koze [22, 23]. Programowanie genetyczne to pró- 
ba zbudowania algorytmu, który będzie tworzył optymal- 
ne z przyjętego punktu widzenia programy (schematy, pro- 
cedury) w oparciu o metodologię algorytmów genetycz- 
nych. Metody programowania genetycznego znalazły na 
przykład zastosowanie przy rozwiązywaniu praktycznych 
problemów z dziedziny tzw. schedulingu, czyli optymal- 


nego zarządzania procesami [1, 8, 19, 26-28, 31], anali- 
zy i klasyfikacji dużych zbiorów danych [35] czy progno- 
zowania zachowania złożonych układów fizycznych [20]. 
Szereg prac jest również poświęconych aspektom teore- 
tycznym stosowania metod programowania genetycznego 
oraz ich łączenia z innymi metodami optymalizacji [4—7, 
10, 11, 17, 29, 30, 32—34]. Osobną grupę stanowią zasto- 
sowania metod programowania genetycznego do rozwią- 
zywania problemów z dziedziny szeroko rozumianej inży- 
nierii naftowej [9, 12]. Szczególne znaczenie mają tutaj za- 
gadnienia związane z optymalną kontrolą i zarządzaniem 
wydobyciem węglowodorów [2, 15, 16]. Problemy te są 
istotne, gdyż decyzja o sposobie eksploatacji złoża ma po- 
ważne konsekwencje finansowe i w dużym stopniu prze- 
sądza o wartości przedsiębiorstw naftowych. 

Tematem prezentowanej pracy jest próba zweryfiko- 
wania możliwości zastosowania algorytmu genetycznego 
(programowania genetycznego) do optymalnego doboru 
wydajności wydobycia ropy naftowej dla prognozowane- 
go czasu wydobycia oraz przyjętego kryterium optymali- 
zacyjnego. Tak postawione zagadnienie optymalizacyjne 
jest przykładem prostego problemu optymalnej kontro- 
li. W artykule przedstawiono podstawy teoretyczne ana- 
lizowanego zagadnienia oraz wyniki działania algorytmu 
optymalizacyjnego realizującego optymalizację dla przy- 
jętego kryterium. 


Podstawy teoretyczne i pojęcia podstawowe 


Niestety nie istnieje żadna spójna teoria systematyzują- 
ca wiedzę na temat algorytmów genetycznych (w tym me- 
tod programowania genetycznego). Do dyspozycji mamy: 

e twierdzenie Hollanda (tzw. twierdzenie o schematach) 
mówiące, ze: 

Wąskie schematy, niskiego rzędu i dobrze dostoso- 

wane rozprzestrzeniają się w kolejnych pokoleniach 

w sposób wyktadniczy; 

e hipotezę „cegiełek” (inaczej „bloków budujących”), we- 
dług której: 

Dobrze dostosowane schematy niskiego rzędu i o ma- 

lej rozpiętości zwane cegiełkami są nieustannie wy- 

bierane, zestawiane i powielane, tworząc ciągi kodo- 

we o potencjalnie wyższym przystosowaniu. 

Powyższe twierdzenia wymagają podania definicji sche- 
matu, jego rzędu oraz rozpiętości. Schemat to w teorii algo- 
rytmów genetycznych wzorzec opisujący podzbiór wszyst- 
kich ciągów podobnych do siebie ze względu na ustalo- 
ne pozycje. Rzędem schematu nazywamy liczbę ustalo- 


nych pozycji w schemacie, natomiast rozpiętość schema- 
tu to odległość między skrajnymi, ustalonymi pozycjami 
w schemacie. 

Podane warunki są zgodne z intuicją. Jeśli schemat jest 
małego rzędu i ma małą rozpiętość, jest mniejsza szansa, że 
w trakcie krzyżowania zostanie rozdzielony. Lepsze przysto- 
sowanie gwarantuje mu natomiast wejście jego reprezentan- 
tów do następnego pokolenia. 

Powyższe twierdzenia tłumaczą mechanizm działania al- 
gorytmów genetycznych jakościowo i nie dają żadnych wska- 
zówek praktycznych. Dla ogólnego przypadku nie istnieją 
dowody ich zbieżności. Twierdzenia takie zostały udowod- 
nione tylko dla prostych funkcji testowych (np. funkcji kwa- 
dratowej) i to w przypadku kodowania binarnego. 

W literaturze poświęconej algorytmom genetycznym (ewo- 
lucyjnym) stosowana jest terminologia mająca swoje źródło 
w teorii ewolucji gatunków. Poniżej zestawiono podstawo- 
we pojęcia stosowane w niniejszej pracy. 

Osobnik — dowolne, możliwe rozwiązanie problemu. 
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Zadaniem algorytmu jest znalezienie w wyniku repro- 
dukcji i rekombinacji (krzyżowania i mutacji) osobnika re- 
alizującego ekstremum globalne. W opisywanym przypad- 
ku osobnikiem jest plik wsadowy do symulatora złożowego 
opisujący w postaci procedury zmiany wydajności wydoby- 
cia ropy naftowej ze złoża. 

Populacja — rodzina osobników. 

Fenotyp — parametry rozwiązania podlegające ocenie. 
Genotyp — kompletny i jednoznaczny zapis parametrów 
osobnika w postaci zakodowanej. 

Gen — pojedynczy element genotypu. 

Chromosom — miejsce przechowywania genotypu osobni- 
ka, czyli jego numeryczna reprezentacja. 

Funkcja przystosowania — funkcja pozwalająca okre- 
Ślić jakość osobnika z punktu widzenia rozwiązywane- 
go problemu. 

Kodowanie rozwiązań — sposób zapisu dowolnego do- 
puszczalnego rozwiązania problemu w postaci genotypu 
osobnika. 

Jakkolwiek w teorii programowania genetycznego znajdu- 
je zastosowanie powyższa terminologia, to jednak ze wzglę- 
du na swoją specyfikę metody te posiadają pewne cechy wy- 
różniające je spośród metod ewolucyjnych: 


e populację stanowią programy komputerowe realizujące 
zadanie optymalizacyjne; 

e funkcja przystosowania ma zadecydować, czy dany pro- 
gram, chociaż w przybliżeniu, wykonuje zadane działa- 
nie. Na przykład, jeśli chcemy stworzyć program, któ- 
ry zawsze będzie znajdywał drogę wyjścia z labiryntu, 
funkcja oceny będzie sprawdzała, jak blisko wyjścia 
program zakończył działanie. Te programy, który znaj- 
dują wyjście, mają najwyższą wartość funkcji przysto- 
sowania, natomiast te, które nie oddalają się od wej- 
ścia, najmniejszą; 

e konieczne jest określenie alfabetu dozwolonych symboli 
i wartości, jakie mogą pojawić się w chromosomie. Mogą 
to być np. liczbowe wartości zmiennych, zmienne, opera- 
tory logiczne, operatory arytmetyczne, słowa kluczowe itd. 
W klasycznym przypadku programowania genetycznego 

programy są przedstawiane jako struktury drzewiaste będą- 

ce graficzną reprezentacją działania programu. Dla naszych 
celów wystarczająca będzie jednak reprezentacja liniowa. 

Główną różnicą pomiędzy klasycznym algorytmem gene- 

tycznym a programowaniem genetycznym jest w tym przy- 

padku fakt, że w tym ostatnim gen może przybierać tylko 
wartości z pewnej grupy wartości dopuszczalnych. 


Sformułowanie problemu 


Cel pracy zrealizowano na przykładzie następującego 
problemu optymalizacyjnego: Poszukujemy takich wartości 
wydajności wydobycia ropy naftowej w przyjętych przedzia- 
łach czasowych, aby w zadanym całkowitym czasie eksplo- 
atacji uzyskać maksymalne wydobycie sumaryczne przy mi- 
nimalnym wydobyciu wody. 

Testy działania zbudowanego algorytmu optymalizacyjne- 
go przeprowadzono dla trójfazowego (ropa, woda, gaz) mo- 
delu złoża PUNQ-S3 (rysunek 1). Model symulacyjny zło- 
ża PUNQ-S3 został udostępniony przez firmę Elf dla celów 
testowych i jest dostępny na zasadach Open Source. Model 
składa się z 2660 bloków (siatka bloków 19 x 28 x 5), z cze- 
go 1761 bloków jest aktywnych. Opisywany niejednorodny 
model symulacyjny charakteryzuje się następującymi wła- 
snościami petrofizycznymi: (1) średnia porowatość — 14,3%, 
(2) średnia przepuszczalność pozioma — 278,8 mD, (3) średnia 
przepuszczalność pionowa — 130,6 mD. W modelu uwzględ- 
niono uskok stanowiący południowo-zachodnią granicę złoża 
oraz dwa aktywne akifery typu Cartera- Tracy ego (przepusz- 
czalność — 137,5 mD, porowatość — 21%, promień — 3000 m, 
miąższość — 19,6 m) okalające złoże od wschodu i południa. 
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Złoże jest eksploatowane trzema odwiertami strefy gazo- 
wej (W-1, W-4, W-12) oraz trzema strefy ropnej (W-5, W-11 
i W-15) zlokalizowanymi wokół kontaktu ropa—gaz. 


Rys. 1. Złoże PUNQ-S3. Widok 3D 
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Schemat działania algorytmu ewolucyjnego 


W celu zbudowania efektywnego algorytmu optymaliza- 
cyjnego bazującego na metodzie programowania genetycz- 
nego konieczne było przyjęcie szeregu założeń dotyczących 
sposobu jego działania. Algorytm zastosowany w pracy opar- 
to na strategii ewolucyjnej zaproponowanej przez Hollanda 
w 1975 r., której schemat zapisany w pseudokodzie przedsta- 
wiono poniżej. W celu uporządkowania opisu poszczególne 
elementy schematu oznaczono liczbami arabskimi. 
Algorytm Hollanda 
begin 
(1) t=0 
(2) inicjalizacja P’ 

(3) ocena P’ 

(4) while [not warunek końca] do 
begin 

(5) T = reprodukcja P“ 

(6) O'= krzyżowanie i mutacja T“ 
(7) ocena O' 


pt = O' 
t=t+l 
end 
end 
gdzie: 


P' — populacja bazowa, 
T' — populacja tymczasowa, 
O' — populacja potomna. 


Algorytm w przedstawionej postaci przetwarza dwie po- 
pulacje: populację bazową P“ oraz populację potomną O'. Wy- 
korzystywana jest również populacja tymczasowa T“ służą- 
ca do przechowywana kopii osobników populacji bazowej. 
W populacjach tych zawarta jest jednakowa liczba osobni- 


ków. W kroku 2 populacja bazowa jest wypełniana losowo 
wygenerowanymi osobnikami. Dla każdego z nich wyzna- 
czana jest w kroku 3 wartość funkcji przystosowania. Po 
przygotowaniu populacji bazowej uruchamiana jest główna 
pętla programu realizującego proces sztucznej ewolucji (4). 
Na początku następuje reprodukcja (5) polegająca na skopio- 
waniu do populacji tymczasowej losowo wybranych osobni- 
ków z populacji bazowej. Prawdopodobieństwa wylosowania 
osobników nie są jednak równe — osobniki o większej warto- 
ści funkcji przystosowania mają większe szanse reprodukcji; 
ponadto stosuje się losowanie ze zwracaniem. W konsekwen- 
cji w populacji tymczasowej w wyniku reprodukcji znajdzie 
się większa liczba kopii lepiej przystosowanych osobników. 
Następnie osobniki z populacji tymczasowej poddawane są 
operacjom genetycznym: krzyżowaniu i mutacji (6). Są one 
kojarzone w pary i dla każdej z nich podejmowana jest de- 
cyzja o krzyżowaniu. Jeżeli jest ona pozytywna, to następuje 
krzyżowanie. W przeciwnym wypadku para nie ulega zmia- 
nie. Następnie każdy z osobników populacji tymczasowej jest 
poddawany mutacji jego genotypu. Decyzja o mutacji genu 
podejmowana jest niezależnie dla każdego genu z zadanym 
arbitralnie prawdopodobieństwem. Utworzone w ten sposób 
osobniki są poddawane ocenie poprzez wyznaczenie warto- 
ści funkcji przystosowania (7) i stanowią populację potom- 
ną, która w kolejnym przebiegu pętli staje się populacją ba- 
zową (4). Pętla jest wykonywana do chwili spełnienia warun- 
ku zatrzymania. Może nim być przeprowadzenie określonej 
liczby iteracji lub znalezienie osobnika o odpowiednio dużej 
wartości funkcji przystosowania (funkcji celu). 

W wersji podstawowej opisany powyżej algorytm po- 
siada dwa parametry: prawdopodobieństwo krzyżowania p, 
oraz prawdopodobieństwo mutacji p,,. 


Kodowanie 


Ponieważ algorytm genetyczny operuje na zakodowanej 
postaci rozwiązania, konieczne jest określenie sposobu jego 
reprezentacji. W ogólnym przypadku metody programowa- 
nia genetycznego do opisu genotypu wykorzystują struktu- 
ry drzewiaste. Dla naszych celów wystarczająca jest repre- 
zentacja liniowa zdefiniowana jako ciąg liczb całkowitych 
identyfikujących wartość wydajności wydobycia dla kolej- 
nych odcinków czasowych na podstawie przyjętego alfabe- 
tu wartości dopuszczalnych. Na przykład chromosom posta- 
ci 1, 3, 4, 2, 5 dla alfabetu wartości dopuszczalnych w posta- 
ci 300, 600, 900, 1200, 1500 oznacza, że złoże eksploatuje 
w kolejnych przedziałach czasowych z wydajnościami 300, 
900, 1200, 600, 1500 m?/dobe. 


Przyjęty sposób kodowania gwarantuje spełnienie po- 
niższych warunków koniecznych do poprawnego działa- 
nia algorytmu: 

e każde rozwiązanie można przedstawić w formie zakodo- 
wanej — w przeciwnym wypadku część rozwiązań była- 
by z definicji nieosiągalna, 

e operacje kodowania i odkodowania są wzajemnie jedno- 
znaczne — w przeciwnej sytuacji możliwe byłoby uzyska- 
nie rozwiązania (genotypu) nieposiadającego interpreta- 
cji fizycznej (fenotypu), 

e operator kodowania nie wprowadza dodatkowych ekstre- 
mów funkcji przystosowania poza tymi, które są odpo- 
wiednikami ekstremów rzeczywistego problemu. 
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Funkcja przystosowania 


Ze względu na przyjęte sformułowanie problemu mak- 
symalizacji w przeprowadzonych testach zbieżności wyko- 
rzystano funkcję przystosowania w postaci: 


BALI 


=o (1) 


gdzie: 

J — współczynnik sczerpania złoża, 
N — sumaryczne wydobycie ropy, 

W — sumaryczne wydobycie wody. 


Reprodukcja 


Operator reprodukcji jest odpowiedzialny za to, które osob- 
niki populacji wejdą do tzw. puli rodzicielskiej, czyli będą pod- 
legać rekombinacji (krzyżowaniu i mutacji), w wyniku czego 
powstanie pokolenie potomne. W niniejszej pracy zastosowano 
reprodukcję rangową, dla której prawdopodobieństwo repro- 
dukcji każdego osobnika jest podane w sposób jawny na pod- 
stawie jego rangi — liczby charakteryzującej jakość osobnika na 
tle populacji. Jest ona sortowana według niemalejącej wartości 
funkcji przystosowania; ranga jest numerem osobnika w takim 


uszeregowaniu. Po określeniu rang definiuje się zmienną loso- 
wą, przypisując każdemu osobnikowi prawdopodobieństwo re- 
produkcji na podstawie jego rangi zgodnie z poniższą formułą: 


1 (Xx) 


Pen Ta @ 


gdzie: 
DP, — prawdopodobieństwo reprodukcji k-tego osobnika, 
r(X,) — ranga k-tego osobnika w populacji. 


Rekombinacja 


Wybór operatorów ewolucyjnych oraz parametrów ich 
działania ma podstawowe znaczenie dla efektywności pra- 
cy całego algorytmu. Trudność polega na tym, że nie istnie- 
ją operatory bezwzględnie ,,lepsze” i „gorsze”. Oznacza to, 
że przyjęty sposób kodowania, który wraz z rozbudowanym 
operatorem selekcji, krzyżowania i mutacji wykazuje bardzo 
dobrą zbieżność dla np. problemu ,,komiwojazera”, może być 
zupełnie bezużyteczny w przypadku kalibracji symulacyjne- 
go modelu złożowego. Specyfika rozwiązywanego problemu 
i wynikająca z niej ograniczona możliwość wykorzystania 
istniejących w literaturze rozwiązań [2, 15], wraz z ilością 
kombinacji możliwych do zastosowania operatorów i war- 
tości parametrów liczbowych, stanowią w tej sytuacji o zło- 
żoności problemu. 

Zastosowane w pracy operatory zostały omówione poniżej. 

Krzyzowanie to proces „mieszania” fragmentów geno- 
typów osobników rodzicielskich w celu wykreowania le- 
piej przystosowanych osobników potomnych z prawdopo- 
dobieństwem p,. W pracy zastosowano metodę krzyżowa- 
nia n-punktowego. Polega ona na wygenerowaniu losowo 
n punktów z zakresu od 1 do liczby genów w chromosomie. 
W efekcie, na podstawie wygenerowanych punktów, każdy 


chromosom jest podzielony na n + 1 fragmentów. Krzyżo- 
wanie polega na budowaniu dwóch osobników potomnych 
przez krzyżowanie dwóch osobników rodzicielskich w wy- 
niku przemiennego kopiowania fragmentów chromosomu 
z osobników rodzicielskich (rysunek 2). 

Mutacja polega na perturbacji genotypu jednego osob- 
nika, przy czym zmiana ta następuje z określonym prawdo- 
podobieństwem p,,. Decyzja o mutacji podejmowana jest dla 
każdego genu każdego osobnika populacji potomnej. Przyj- 
muje się ponadto, że niewielkie perturbacje są bardziej praw- 
dopodobne niż duże. 


Rys. 2. Zasada krzyżowania n-punktowego 


Przykłady optymalizacji 


Adaptacja parametrów 
Jednym z głównych problemów wynikających ze stoso- 
wania algorytmów genetycznych, a w szczególności metod 
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programowania genetycznego, jest wystepowanie ,,przed- 
wezesnej zbieznosci” (bedacej niczym innym jak szybkim 
osiągnięciem przez algorytm pewnego ekstremum) i ,,zablo- 


kowania” dalszego procesu optymalizacji. W odróżnieniu od 
metod klasycznych algorytmy genetyczne posiadają mecha- 
nizmy zapobiegania tego rodzaju sytuacjom. Prawie zawsze 
stosowana jest mutacja umożliwiająca wprowadzenie loso- 
wego, ale kontrolowanego zaburzenia procesu krzyżowa- 
nia. Drugim wykorzystywanym często zabiegiem jest ada- 
ptacja prawdopodobieństw krzyżowania p. i mutacji p,. Po- 
lega ona na tym, że w każdym pokoleniu p. i p, są wyzna- 
czane od nowa, w zależności od wartości wybranego para- 
metru charakteryzującego populację. 

W prezentowanym algorytmie zastosowano adaptację ste- 
rowaną współczynnikiem zróżnicowania Z,;: 


fmin(X) 


Ta = can ©) ©) 


W powyższej formule przez f(X) oraz f ean(X) oznaczo- 
no odpowiednio najgorsze i średnie przystosowanie w popu- 
lacji. Parametr ten jest tym większy, im większe jest zróżni- 
cowanie osobników. W granicznym przypadku, gdy popula- 
cja składa się z identycznych osobników, Z = 0. 


artykuły 


Dla każdego pokolenia prawdopodobieństwa krzyżowania 
p.i mutacji p, wyznaczane były według następujących formuł: 


Pe = Pe min + Za 3 (p. max Pe ran) (4) 


Pm = Pm max Za j (Pm max —Pm an) (5) 


gdzie indeksy max, min odnoszą się do zadanych maksymal- 

nych i minimalnych wartości parametrów. 

Działanie tak zdefiniowanego mechanizmu adaptacji jest 
następujące: 

e jeżeli zróżnicowanie populacji rośnie (Z —> 1), automa- 
tycznie zwiększane jest prawdopodobieństwo krzyżowa- 
nia, a zmniejszane prawdopodobieństwo mutacji, 

e jeżeli zróżnicowanie populacji zmniejsza się (Z — 0), 
mechanizm działa odwrotnie, tzn. zmniejsza prawdopo- 
dobieństwo krzyżowania, a zwiększa prawdopodobień- 
stwo mutacji. 

Przedstawione rozwiązanie z jednej strony gwarantuje 
różnorodność populacji, z drugiej zaś zapewnia dobrą eks- 
ploatację przestrzeni rozwiązań. 


Wyniki symulacji 


Zbudowany algorytm ewolucyjny zaimplementowano 
w kompilatorze Parallel Studio XE 2011 firmy Intel. Przy 
wyznaczaniu wartości funkcji przystosowania użyto symu- 
latora złożowego ECLIPSE 100 firmy Schlumberger. 

Dla każdego z przedstawionych poniżej przykładów opty- 
malizacji przyjęto: 
e maksymalna liczba generacji: 35, 
e liczba osobników w pokoleniu bazowym: 10. 

Ponadto dla prawdopodobieństwa krzyżowania p. oraz 
mutacji p, przyjęto następujące wartości graniczne: 


Pe min = 0,5; Pe max iw 0,9 
Pm min = 0,0 1 ; Pm max = 0,3 


oraz formuły opisujące ich zmienność (3—5). 

W trakcie przeprowadzonych testów wykonywano 30-let- 
nie prognozy wydobycia, maksymalizując sumaryczne wy- 
dobycie ropy naftowej ze złoża zgodnie z formułą (1). 

Testy prowadzono dla trzech różnych długości chromo- 
somu dla krzyżowania |-punktowego (test — 1): 

e długość chromosomu: 5 (5 przedziałów czasowych po 

6 lat każdy), 

e długość chromosomu: 10 (10 przedziałów czasowych po 

3 lata każdy), 

e długość chromosomu: 15 (15 przedziałów czasowych po 

2 lata każdy), 
oraz dla 1, 2 13 punktów krzyżowania dla długości chromo- 
somu równej 10 (test 2). 


Zbiór dozwolonych wydajności wydobycia ze złoża skła- 
dał się z 6 wybranych arbitralnie wartości: 300, 600, 900, 
1200, 1500, 1800 m*/dobe. 

Zadaniem zbudowanego algorytmu był taki dobór wy- 
dajności wydobycia w poszczególnych przedziałach czaso- 
wych, aby uzyskać maksymalne wydobycie ropy przy mi- 
nimalnym wydobyciu wody. Na podstawie wyników prze- 
prowadzonych symulacji przeanalizowano wybrane aspek- 
ty jego działania. 

Rysunek 3 przedstawia zależność wyników optymaliza- 
cji od długości chromosomu dla 3 wybranych jej wartości 
(5, 10, 15). W praktyce oznaczało to podziały 30-letniego 
okresu prognoz na 6-, 3- i 2-letnie okresy, w których wydaj- 
ność wydobycia była stała. Testy w tym przypadku prowa- 
dzono dla krzyżowania |-punktowego. Dla każdego z przed- 
stawionych przypadków uzyskano identyczną wartość funk- 
cji przystosowania, przy czym nie stwierdzono jednoznacz- 
nej zależności szybkości procesu optymalizacji od długości 
chromosomu. Charakterystyczne jest to, że spośród trzech 
długości chromosomu najgorszą zbieżność otrzymano dla 
najkrótszego chromosomu odpowiadającego najdłuższym, 
czyli 6-letnim przedziałom czasowym. 

Na rysunku 4 przedstawiono zależność zbieżności algo- 
rytmu od liczby punktów krzyżowania dla krzyżowania 1-, 
2- i 3-punktowego (test 2). Podobnie jak dla testu zależności 
zbieżności od długości chromosomu (rysunek 3) nie stwier- 
dzono zależności od badanego parametru. W każdym z przy- 
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Rys. 3. Zależność zbieżności algorytmu od długości 
chromosomu 
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Rys. 4. Zależność zbieżności algorytmu od liczby punktów 
krzyżowania 


padków (1, 2, 3 punkty krzyżowania) uzyskano identyczne 
maksimum, zgodne z wynikiem testu 1. 

Należy zauważyć, że brak widocznej zależności pomię- 
dzy zbieżnością a badanymi parametrami (długość chromo- 
somu, liczba punktów krzyżowania) wynika ze stochastycz- 
nego charakteru zaproponowanej metody optymalizacji. 


Weryfikacja takiej zależności wymaga powtórzenia procesu 
optymalizacji dla każdej wartości badanego parametru kilka 
tysięcy razy i porównania uśrednionych wyników. 

Rysunek 5 przedstawia sumaryczne wydobycie płynów 
złożowych dla rozwiązania najkorzystniejszego ze wzglę- 
du na kryterium optymalizacyjne (1). Dla najlepszego roz- 
wiązania otrzymano stałą wydajność wydobycia w całym 
okresie prognozy na poziomie najniższym spośród dozwo- 
lonych — 300 m*/dobe. Rozwiązanie na najniższym możli- 
wym poziomie wydajności wynika z postaci funkcji przysto- 
sowania, która uwzględnia zarówno wydobycie ropy nafto- 
wej, jak i wody złożowej. Fakt, że w miarę wzrostu wydaj- 
ności wydobycia następuje gwałtowny, nieliniowy przyrost 
wydobycia wody, powoduje, że dla zadanej postaci funkcji 
przystosowania algorytm promuje jak najniższe wydajności 
wydobycia, dla których woda pojawia się późno i w małej 
objętości. Liniowa charakterystyka sumarycznego wydoby- 
cia ropy naftowej wynika z faktu, że dla optymalnego roz- 
wiązania w każdym przedziale czasowym otrzymano takie 
same wartości wydajności wydobycia ropy naftowej. 
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Rys. 5. Sumaryczne wydobycie płynów złożowych 


Podsumowanie 


Zadaniem artykułu było zweryfikowanie możliwości za- 
stosowania metod programowania genetycznego dla celów 
optymalizacji wydobycia ropy naftowej przy założeniu mi- 
nimalizacji wydobycia wody złożowej. W wyniku wykona- 
nych symulacji potwierdzono możliwość zastosowania wy- 
branej metody do optymalizacji wydobycia węglowodorów 
przy założonych warunkach ograniczających. 

Przeprowadzone analizy symulacyjne pozwalają sformu- 
łować następujące wnioski i uwagi: 

1. Zastosowana metoda okazała się być efektywna dla roz- 
wiązania postawionego problemu optymalizacyjnego. 
Przeprowadzone testy wykazały zadowalającą zbieżność 
zarówno dla różnych długości chromosomu, jak i róż- 
nej liczby punktów krzyżowania. Co więcej dla każdego 
z przeprowadzonych testów otrzymano w każdym z ana- 


766 Nafta-Gaz, nr 10/2017 


lizowanych wariantów identyczne rozwiązanie, co świad- 
czy o bardzo dobrej powtarzalności działania algorytmu. 

2. Przyjęcie funkcji przystosowania w postaci 1 powoduje pre- 
ferowanie przez algorytm optymalizacyjny niskich wydaj- 
ności wydobycia ropy ze względu na szybki przyrost wy- 
dobycia wody w przypadku zwiększania wydobycia ropy. 

3. Operator mutacji ma kluczowe znaczenie dla efektyw- 
nego przeszukiwania przestrzeni rozwiązań. Przeprowa- 
dzone symulacje wskazują, że zastosowanie prawdopo- 
dobieństwa mutacji na zbyt wysokim poziomie powodu- 
je problemy z utrzymanie zbieżności algorytmu. 

4. Wskazane jest stosowanie adaptacji parametrów sterują- 
cych działaniem algorytmu. Wynika to z potrzeby różni- 
cowania osobników w zależności od stanu procesu opty- 
malizacji. 
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